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The relation between the geometric phase and quantum phase transition has been discussed in the 
Lipkin-Meshkov-Glick model. Our calculation shows the ability of geometric phase of the ground 
state to mark quantum phase transition in this model. The possibility of the geometric phase or its 
derivatives as the universal order parameter of characterizing quantum phase transitions has been 
also discussed. 



PACS numbers: 75.10.Pq, 03.65.Vf, 05.30.Pr, 42.50.Vk 



O 
O 

(N 

s ■ 
< 



> 

m 
O 

oo 

o 

O 

S3 

^— > ■ 
G ■ 
cd 

3 : 



I. INTRODUCTION 

Recently the understanding of quantum phase tran- 
sition has emerged from the fundamentals of quan- 
tum mechanics, especially from the entanglement point 
of view [31 ■ The ground state structure have been shown 
to be affected appreciably by the critical points, as illus- 
trated initially in the one-dimensional Ising model with 
transverse magnetic field ■ Intriguing by these pioneer- 
ing works, agreat deal of efforts have attributed to this 
region [1,01|(J. It originates from the belief that quantum 
phase transition should be connected to the entanglement 
or its derivatives in many-body systems Q. However 
the correspondence between the entanglement and quan- 
tum phase transition in many-body systems is ambiguous 
0, , that is because of the absence of the proper mea- 
surement of entanglement in many-body systems. On the 
other hand, geometric phase [9j as a measurement of the 
curvature of the Hubert space, has been first connected 
with the quantum phase transition in the one-dimension 
spin-chain system |lfjj , in which the topological character 
of the relative geometric phase between the first excited 
state and the ground state have been shown the ability 
of detecting the critical points. Furthermore, Zhu ex- 
tended this study under the thermodynamic limit and 
found that the geometric phase of the ground state in 
the one-dimensional XY model was non-analytical and 
its derivative with the coupling constant was divergent 
closed to critical points [Tjj . Moreover the scaling behav- 
ior of the geometric phase of ground state clearly distin- 
guished two different types of quantum phase transitions 
at the critical point in the one-dimensional XY model. 
Consequently the general theory about the relation be- 
tween geometric phase of the ground state and the critical 
point has been constructed |12| , in which the topological 
property of Berry's loops including critical points in the 
parameter space has been shown the ability of detecting 
the critical points. 

Although the theory is rapidly developing, the verifi- 
cations of examples only focus on the one-dimension XY 
spin-chain systems. The reason is that XY model, which 
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have been shown first-order phase transition closed to the 
critical point, can be converted into spinless fermionic 
system by Jordan- Wigner transformation, and the en- 
ergy spectrum can be determined exactly 13] . However 
the interacting between spins in this model is short-range 
(the neighbor-nearest coupling) and the anisotropy of 
Heisenberg interaction is indispensable for the construc- 
tion of long-range order Recently the Berry phase 
in Dick model has been examined under thermodynam- 
ics limit ^4|. The authors showed that the Berry phase 
displayed the non-analyticity close to critical point where 
this model exhibits a second-order phase transition, and 
its derivative with the coupling constant was also dis- 
continued and showed a cusp closed to this point. The 
discontinuity of Berry phase is very similar to that in 
XY model, but we should point out that the two model 
belong to different orders of quantum phase transitions 
respectively. In this paper, we will show a special situ- 
ation, in which the geometric phase of the ground state 
behaves differently and itself distinguishes the different 
quantum phase transitions. 

In particular we study the geometric phase in a sys- 
tem of spins with a collective coupling described by the 
Lipkin-Meshkov-Glick (LMG) model 1 1 5| . which first in- 
troduced forty years ago in nuclear physics. This model 
has received much attention because of his apparent sim- 
plicity and popularity in the past. It provides a simple 
description of the tunneling of bosons between two de- 
generate levels and can thus be used to describe many 
physical systems such as two-mode Bose-Einstein con- 
densates or Josephson junctions ^?|]- More recently 
the entanglement in this model has received great atten- 
tion because of available numerical calculations and plen- 
tiful phase diagram || 0] . Under the thermodynamics 
limit, its phase diagram can be simply established by a 
semiclassical approach |19| . For large particle number 
N but finite, the situation is complicated and the nu- 
merical analysis was implemented using the continuous 
unitary transformations |2Cj . The significant difference 
between this model and XY model is the long-range in- 
teraction, and the system cannot be converted into the 
spinless fermionic system. Hence it is of great interest to 
study the relation of geometric phase and phase transi- 
tion in this model. As will display in the remaining of this 
paper, the geometric phase of ground state in this model 
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behaves differently and reflects faithfully the existence of 
the critical points. 

The paper is organized as follows. In Sec. II, we first 
introduce the LMG model, and in the limit of large N 
but finite, we obtain the ground state analytically by the 
Holstein-Primakoff representation and calculate the ge- 
ometric phase. In order to show the university of our 
results, in Sec. Ill, we study the phase transition in a 
more complex situation. It is surprising that the geo- 
metric phase of ground state detects very rigorously the 
critical points. Finally we discuss the implications of our 
results and the differences from the XY model. 



II. THE LIPKIN-MESHKOV-GLICK MODEL: 
BIAXIAL CASE 

The LMG model describes a set of N spins half coupled 
to all others with a strength independent of the position 
and the nature of the elements and a magnetic field in 
the z direction. The Hamiltonian can be written 



H = -±(S 2 x+1 S 2 y )-hS z , 



(1) 



in which S a = J2i=i a a/^( a = x i z ) an d the a a is the 
Pauli operator, N is the total particle number in this 
system. The prefactor 1/N is essential to ensure the 
convergence of the free energy per spin in the thermo- 
dynamic limit. For any anisotropy parameter 7 € [0, 1], 
the Hamiltonian Q preserves the total spin and does 
not couple the state having spin pointing in the direction 
perpendicular to the field, namely 
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[tf,S 2 ]=0, [H,l[ai]=Q. 



(2) 



An important character that differentiates LMG model 
from the XY model is the long-rang interaction between 
particles, which induces a second-order phase transition 
at h = 1 when N — > 00 |19j. 

The diagonalization of Eq. can be obtained by 
introducing the Holstein-Primakoff representation of the 
spin operator and then truncate the resulting bosonic 
Hamiltonian to lowest order 20]. Consequently we diag- 
onalize it thanks to the Bogoliubov transformation. The 
first thing is to perform a rotation of the spin operators 
around the y direction, that makes the z axis along the 
so-called semiclassical magnetization [2(j in which the 
Hamiltonian Eq. has the minimal value in the semi- 
classical approximation. This can be done as 




cos6» -sin<9 \ / S x 
1 ] ( S„ 
sin 9 cos 9 



(3) 



in which 9 — for h > 1, 9 — arccos/i for < h < 1 
Then the Hamiltonian Eq. is converted into 

rr 1 r • 2 fl( J 2 , COS 2 + 7 ~ 2 ~2 

H = --[sin 9S Z + (S -S z )] 

-{b^b z + b bz+h.c.j — 



in which = S x ± iS y . 

In order to obtain the geometric phase of ground state, 
we consider the system has a rotation g(4>) around the 
new z direction. The Hamiltonian becomes 



N 

h z cos 9S Z - sm6{S + + S-), 



AN 



H{<j>) = ~g{4>)Htf{<i>) 



(5) 



in which g{4>) = e lc ^ Sz . Then we can use the Holstein- 
Primakoff representation, 



S + (<p) = {N-a^a) 



(6) 



in which a^) is bosonic operator. Since the z axis is 
along the semiclassical magnetization, a^a/N < 1 is a 
reasonable assumption under low-energy approximation, 
in which N is large but finite. Keeping the terms of order 
A^iV 1 / 2 ,^ , Eq. © becomes 

H{4>) =Ne + Aa^a + r(a t2 e- 2 '^ + a 2 e 2 ^), (7) 

in which 

1 



e = — -(sin 2 9 + 2/icos6>), 

a • 2 n 7 + cos 2 a 
A = sm 9 - 



h cos t 



r 



7 — cos 



(8) 



Obviously the equation above can be diagonalized by the 
standard Bogoliubov transformation, 

b{4>) = cosha;ae # + sinhiraV 40 (9) 

Then the Hamiltonian is 

H diag {<t>) =Ne + <r + ^b^^bCP), (10) 

in which, 




2-/i 2 -7' 



< h< 1 



It should be careful about the physical interpretation of 
A D , which may not describe true gap of the system [2(J . 

Now it is time to determine the ground state | <?(</>)), 
which can be obtained by applying the relation, 



&(</>) \g{<t>)) =0 



(12) 



Substituting Eq. © into the equation above, one obtains 
the ground state, 



(4) 
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(13) 
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FIG. 2: <p s [Arc] vs. ft with different particle number N. 
We have chosen 7 = 0.5 (a) and 7 = 1 (b) for this plot. The 
dashed and solid lines correspond respectively to N = 4, 1000. 



FIG. 1: The geometric phase ip g [Arc] vs. the anisotropic 
parameter 7 and ft. For specification, we have chosen the 
summation from to 100 in the expression of ip g (Eq. 1141 '). 
The divergent character of tp g is clearly displayed at ft. — > 1 in 
this figure. 



in which n!! = n(n — 2)(n — 4) • • • and rill = 1 
for n < 0. \n) is the Fock state of bosonic op- 
erator a^) and the normalized constant is C 2 = 
V^/ 2] 2 sinh 2 x^f^ tanh 2 ^" 1 ) x. One should note 
that in order that the summation is convergent, 
I tanha;| < 1. 

The geometric phase of the ground state, accumu- 
lated by changing <f> from to tt, is determined by 
(p g = —ij^d(j>(g((f>)\d < p\g((t>)). The direct calculation 
shows 



(fig = 7T[1 



El^n^tanh 2 ^). 
E^^tanh 2 ^ 1 ), 



(14) 



which have been plotted with 7 and h in Fig. ^ It is 
obvious that tp g is divergent at the point h — 1 where 
the LMG model has been proved to experience a second- 
order phase transition, independent of the anisotropy 
|l9j - This divergency has never been explored in the 
previous studies [l(J, HH Il2> UH an d shows distinguished 
character from the XY and Dicke models. In fact, since 
the adiabatic condition has been destroyed when h — > 1 
because of the degeneracy between the ground state and 
the excited state and the degeneracy enhances geomet- 
ric phase 0, this phenomenon is natural. However, the 
essential reason is the collective interaction in the LMG 
model, which is absent in the XY model and makes the 
long-range correlations in this system. 

The finite size effect is also examined in this case by 
choosing different N, which can be shown in the Fig. |2Ja) 
and (b). In this figure, we only draw for N = 4, 1000 
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FIG. 3: The scaling behavior of tp g [Arc] vs. iV with 7 = 0.5 
when h — > 1. The slope of the line is close to —1. 



and for higher values of N, the curves almost do not 
change. From the figures we note that the divergency of 
(fig seems to be insensitive to N. However, since phase 
transition only happens under thermodynamic limit, this 
phenomenon attributes to the ground-state degeneracy 
at the critical point. 

The scaling behavior of <p g plotted with N has also 
been explored in Fig. [21 It is obvious that one can obtain 



-N. 



(15) 



Furthermore the scaling is independent of 7, which means 
that for different 7, the phase transitions belong to the 
same university class. We also check the scaling behav- 
ior when 7 = 1 and find that this case has the same 
behavior as displayed in Fig. |3| It means that the phase 
transition for 7 = 1 is the same university class as that for 
7^1. This phenomenon is different from the XY model, 
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in which the isotropic and anisotropic interactions be- 
long respectively to different university classes [Tlj . that 
comes from the collective interaction in the LMG model. 



III. UNIAXIAL MODEL 

In this part, we will discuss the phase transition in a 
more complex system, which is the generalization of the 
LMG model Eq. The Hamiltonian can be written as 



H — — 77 S~ 

N x 



hxS x h z S z 



(16) 



with h z > 0. The phase diagram of this model is ob- 
viously dependent on the both parameters h x ,h z , and 
moreover a proper order parameter characterizing the 
transition is difficult to build. Recently the correspon- 
dence between this model and a two-level boson prob- 
lem introduced in nuclear physics has been constructed, 
which permits one to get an order parameter plj . The 
phase transition is then clear in this model; a first-order 
transition occurs at h x — with h z < 1 and a second- 
order one occurs at h x — with h z = 1. For h z > 1 or 
h x ^ 0, no transition is found. 

Now we will try to determine the phase transition by 
calculating the geometric phase of the ground state. Sim- 
ilar procedures as in previous section can be applied 
for this purpose. The first step is to make the sys- 
tem have a rotation around the z direction and then 
the Hamiltonian 1161 becomes H(<f>) — g(4>)Hg* (</>). Next 
step is to introduce the Holstein-Primakoff transforma- 
tion Eq. JBJ). One should note that the approximation 
a^a/N <C 1 is invalid in this case since we cannot find the 
semiclassical magnetization. However a simple canoni- 
cal transformation can be used to resolve this problem, 
a (t) e (-W = &(t) e (-W + y/NX with |A| < 1. This trans- 
formation provides a macroscopic expectation value of S z 
which is order of N, and then one has b^b/N -C 1. Under 
the limit that N is large but finite, it is enough to expand 
H((j)) to the order N°. After the exhausted calculation, 
H(cj)) is written as, 

H(<f>) = e (A) + n(A)(6e** + tie'*) + 

r(A)(6 2 e 2 ^ + 6 t2 e~ 2 ^) + A(A)6 f 6, (17) 

in which, 



e (A) = -iV[^(l-2A 2 ) + A 2 (l-A 2 ) + Ml-A 2 )] 

2-1 
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r(A) 
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FIG. 4: The geometric phase ip g [Arc] vs. h x and h z . We 
have chosen the summation from to 100 in the expression 
of ipg Eq. dJ, and for the convenience of viewport, we have 
drawn for — ip g in this plot. 



The crucial step is to choose Ao properly in order that the 
linear term (the second term in Eq. (|17fl ) is vanishing. 
It can be realized by solving the following equation, 



X h z 



h x (l~2X 2 ) 



2 VI 



A 2 , 



-A (1-2A 2 ) = 0. 



(19) 



'4(l_A2)3/2 



(18) 



In fact the equation above can be reduced into the bi- 
quadratic equation (h z — y) 2 (l — y 2 ) — h x y 2 = with 
y = 1 — 2Aq, which can be solved numerically. 

Substitute Ao into Eq. i|17|) an d then one get the 
quadratic Hamiltonian, 

H = e (Ao) +r(A )(6 2 e 2 ^ + 6 t2 e~ 2 ^) + A(A )6 t 6, (20) 

which obviously could be diagonalized by the standard 
Bogoliubov transformation. Consequently the geometric 
phase of ground state can be determined directly, which 
has the same form as Eq. (|14l) , but different definition of 

x, determined by the equation tanh2x = ^t^~j ■ 

A schematic demonstration of the geometric phase is 
presented in Fig. It is obvious that there are two 
regions divided by h z — 1, and the geometric phase is 
divergent at point h x = 0, h z = 1. In the region h z < 1, 
geometric phase is non-analytical at point h x = 0, which 
means the appearance of phase transition, and in the 
other region, the geometric phase is the smooth function 
of h z , h x and no phase transition is found from the fig- 
ure. This phenomenon is consistent with the conclusion 
in Ref. but in our calculation we do not need to 
find a proper order parameter to characterize the phase 
transition and the geometric phase of ground state faith- 
fully marks these transitions. A detailed demonstration 
is also provided in Fig. [3J From the figures, one easily 
finds that the geometric phase tp g has a cusp at h x = 
for h z < 1 (see Fig. Ufa)), which implies the first-order 
phase transition. Furthermore ip g is divergent at h. x — 
when h z — 1 (see Fig. Efb)), which is similar to that in 
the standard LMG model (see Figs. ^,nd|21) an d means 
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FIG. 5: The geometric phase <p g [Arc] vs. h x with different 
h z . We have chosen h z = 0.5 (a), h z = 1 (b) and h z = 2 (c) 
for this plot. The figure (b) has been compressed because of 
the divergency of the value tp. 

that there is a second-order transition. For h z > 1 (p g is 
the smooth function of h x , h z and there is no transition 
(see Fig. GJc)). 

It is surprising that the geometric phase of ground state 
itself can differentiate phase transition without need of 
introducing a proper order parameter, which may be dif- 
ficult to find. Thus it is a natural speculation that the 
geometric phase of ground state can serve as a universal 
order parameter. Further discussions will be presented 
in the next part. 

IV. CONCLUSIONS AND DISCUSSIONS 

The geometric phases of ground state in the LMG 
model and its generalization have been discussed in this 
paper. Our calculations show that the geometric phase 
faithfully reflects the phase transition in this model; when 
there is a second-order phase transition, the geometric 
phase behaves divergent (see Figs. nand|S{b)). However 
when there is a first-order transition, non-analyticity of 
geometric phase appears at the critical point (see Fig. 
Eta)). These phenomena may come from the degeneracy 
of the ground state in the system, which usually induces 
the mixture of different phases. The different behaviors 



of geometric phase closed to critical points could origi- 
nate from the different excitations (i.e. gapped or gap- 
less). Furthermore we find that the energy of ground 
state is not a good parameter of marking the phase tran- 
sition, because the energy of ground state is degenerate 
independent of the phase transition is second-order or 
first-order. Recently Tian and Lin have shown that the 
continuous quantum phase transition are actually caused 
by level crossing of the low-lying excited states of the sys- 
tem [2^ • This conclusion also shows that the only energy 
of the ground state is not suitable for the characterization 
of quantum phase transition. 

This leads to a question what physical quantity is suit- 
able for characterizing the qua ntum phase transition. 
With respect of the work Ref and our calculations, it 
seems to provide us an hint that the geometric phase of 
ground state or its derivatives could serve as an univer- 
sal order parameter to characterize different phase tran- 
sitions. This speculation is natural since the geomet- 
ric phase faithfully measure the curvature of the Hilbcrt 
space (or phase space for classical mechanics) and the 
broken of symmetry of Hilbert space must be reflected in 
the geometric phase. 

Another aspect of importance is the differences be- 
tween our model and the one-dimensional XY model. 
The crucial point is the collective interaction in LMG 
model, which is absent in the XY model. A main re- 
sult of this interaction is that the LMG model cannot 
be converted into the spinless fermion system. This fact 
makes the conclusion different from that in Ref. 01 > m 
which the topological behavior of Berry phase can detect 
the critical point. However, in this model, the geometric 
phase is divergent when there is second-order transition 
and the detection maybe invalid. 

In conclusion the relation between the geometric phase 
of ground state and the phase transition in LMG model 
has been constructed in this paper. Different from the 
results in the one-dimensional XY model, the singularity 
of the geometric phase itself can serve as the signature 
of critical points. Moreover we discuss the possibility 
of the geometric phase of ground state or its derivatives 
serving as the universal order parameter to characterize 
the quantum phase transitions. 

This work was supported by NSF of China under 
grants 10305002 and 60578014. 
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